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The attosecond, time-resolved X-ray double-quantum-coherence (XDQC) four wave mixing sig- 
nals of formamide at the nitrogen and oxygen K-edges are simulated using restricted excitation 
window time-dependent density functional theory (REW-TDDFT) and the excited core hole (XCH) 
approximation. These signals, induced by core exciton coupling, are particularly sensitive to the 
level of treatment of electron correlation, thus providing direct experimental signatures of electron 
core hole many-body effects and a test of electronic structure theories. 



I. INTRODUCTION 

New attosecond X-ray table top sources [H, [If and the 
X-ray free electron lasers (XFEL) [sl, 0] allow to extend 
nonlinear spectroscopic techniques, originally developed 
to study valence excitonic systems, to molecular core ex- 
citations. XFEL sources are bright enough to saturate 
the core-excitation transitions, [5, 6] and should enable 
fundamental questions about the quasiparticle descrip- 
tion of many-electron states to be experimentally ad- 
dressed. Nonlinear techniques designed to probe quan- 
tum correlations in systems at lower energies with slower 
dynamics may thus be extended to the X-ray frequency 
regime. The core electrons are tightly bound to specific 
atoms and are strongly coupled to a slower bath of cor- 
related valence electrons shared between atoms. ID- and 
2D- stimulated X-ray Raman techniques have been pro- 
posed to examine valence electron dynamics using core- 
holes as ultrafast switches. 

This paper focuses on doubly core-excited states 
(DCESs), in which the main players are the static core 
holes and the virtual orbitals. These states correspond 
to the high-energy limit of double valence excitations, 
which are challenging to treat with density functional 
theory (DFT) methods. [7] Multiply excited states can 
be prepared by short and intense laser pulses. [8] It has 
been suggested that DCESs carry information about 
the chemical environment of a selected atom in double 
photoionization core-excitation spectroscopv. j9l-[ll| Re- 
cent experimental [12., J^3] and theoretical [9, fll| work 
demonstrated a static frequency shift of doubly-core- 
photoionized states that depend on their local chemi- 
cal environment. The double-quantum-coherence (DQC) 
signal preferentially targets doubly excited states in 
an excitonic system. [15| This signal strongly depends on 
the coupling between excitons; infrared DQC signals, 
which have been used to study vibrational couplings in 
the amide I band,[lg] sensitively depend on anharmonic- 
ities, vanishing for harmonic vibrational systems. In the 
visible, DQC has been used to investigate the excitations 
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of a dye molecule in ethanol,[T75 and theoretical stud- 
ies have demonstrated a strong dependence on electron 
correlations in model two and three-level semiconductor 
systems [isj. It can monitor the breakdown of mean- 
field theories of electron correlation in molecular systems. 
This technique has also been applied to measure spin se- 
lected biexciton coupling using circularly polarized pulses 
[2Q|, and interatomic couplings in a dilute potassium 
vapor [215. The X-ray version of this technique (XDQC) is 
sensitive to correlation and exciton scattering in DCESs, 
making it an attractive experimental test for many-body 
electron structure techniques for strongly correlated sys- 
tems. [22-24] XDQC detects only systems in which the 
singly core excitations affect each other and the doubly 
core-excited wavefunction may not be factorized into an 
outer product of two singly core-excited wavefunctions. 
For uncoupled core excitations, two contributions to the 
signal (see Sec. lIVp with opposite signs cancel out and 
the signal vanishes. Incomplete cancellation therefore 
provides a sensitive measure of correlation between core 
excitations. 

Core-excited states can be calculated at various lev- 
els of theory. The XDQC signal was first simulated 
for all nitrogen excitations of isomers of aminophenol 
with a simple equivalent core approximation (ECA) de- 
scribes the effect of the core hole on the valence ex- 
cited states. [25] In this study we employ restricted ex- 
citation window time-dependent density functional the- 
ory (REW-TDDFT), [lilai] a response formalism that 
incorporates the core hole valence coupling through 
an exchange-correlation functional. This approach was 
recently used to calculate the single core excitations 
and two dimensional stimulated X-ray Raman signal of 
cysteine. |33| 

Here we simulate the XDQC signals of formamide, 
a small organic molecule containing carbon, nitrogen 
and oxygen, which is often used as a model for the 
peptide backbone in proteins (see Fig. [3|). Its X-ray 
spectra has been investigated both experimentally [3J- 
|36| and theoretically. [37] Formamide has been used as 
a benchmark in theoretical studies of double core hole 
spectroscopy [m |39| and in molecular dynamics [40| and 
vibrational spectroscopy applications. |4lll42j 

We summarize the theoretical methods and computa- 
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tional challenges of double excitations in Sec. [Ill Our 
treatment of core and valence excited states is presented 
in Sec. HITl Sum over states XDQC expressions for sig- 
nals involving the core-excited eigenstate frequencies and 
transition dipole matrix elements are given in Sec. IIVI 
We then apply the REW-TDDFT to calculate the ori- 
entationally averaged XANES and XDQC signals of for- 
mamide at the N and O K-edges. The results are dis- 
cussed in Sec. |Vll 



II. ELECTRONIC STRUCTURE SIMULATIONS 
OF DOUBLE EXCITATIONS 

Electronic excited states can often be represented as 
linear combinations of single excitations, where a sin- 
gle electron is promoted from an occupied to a vir- 
tual orbital. Some of the most popular and inexpensive 
computational methods for excited states, such as con- 
figuration interaction singles (CIS) and time-dependent 
density functional theory (TDDFT), are based on this 
picture. While often adequate for molecules or high 
bandgap materials they fail badly for double excita- 
tions in conjugated or metallic systems with strong va- 
lence band correlations. Double-exciton states, which 
are common in molecular crystals (43j and materials with 
strong spatial localization, [4^ also play a role in coni- 
cal intersections, long-range charge-transfer excitations 
and autoionizing resonances. Low- lying excited states 
in polyenes with significant double excitation charac- 
ter are notable examples in which the single excitation 
picture fails qualitatively to describe the svstem. (45l452| 
Strong double-excitation features in X-ray absorption 
near edge structure (XANES) spectroscopy have been 
reported for ferrocene and ferrocenium compounds, [53| 
and it has been shown that double or higher order 
excited configurations are necessary to construct the 
spin-symmetry-adapted wave- functions of molecules with 
open-shell ground states. [54- 57] 

Formally any state with energy close to the sum of 
two single excitation energies can be considered a dou- 
ble excitation. However this may not be always justified. 
The doubly excited states in polyenes mentioned above 
are examples where this assignment completely breaks 
down, since they are lower in energy than any single ex- 
cited state. [48] To define double excitations, we must first 
specify the reference single-particle theory and define the 
orbitals and their energies. In some cases a Hartree-Fock 
ground state reference suggests a significant double exci- 
tation character for an excited state while a Kohn-Sham 
ground state reference does not.0, Various high-level 
ab initio methods also show very different double excita- 
tion character, [48, 58] since DCESs are directly related 
to correlation. The XDQC technique creates DCESs by 
using two pulses. In contrast with valence shake-up ex- 
citations in XANES induced by the core hole Coulomb 
interaction, in two photon core excitation each core is 
prepared by a different photon. In this paper we consider 



the valence relaxation in the field of the first core hole 
explicitly, and neglect any accompanying Auger, nuclear 
or electronic processes within the very short duration of 
the measurement. 

High-level ab initio techniques are usually required to 
handle double excitations in valence band excitations. (59| 
These include complete active space self-consistent 
field (CASSCF), complete active space perturbation 
theory of second order (CASPT2),[5i| coupled clus- 
ter (CC),|60] multireference configuration interaction 
(MRCI),[52, 61] symmetry-adapted cluster configura- 
tion interaction (SAC-CI),(62| algebraic diagrammatic 
construction (ADC)[48'-f50j and multireference M0ller- 
Plesset perturbation theory (MRMP)[46, 58]. These ac- 
curate techniques are computationally expensive and lim- 
ited to small molecular systems such as polyenes with 
several carbon atoms. 

TDDFT [^,[63 balances accuracy and computational 
cost for excited states. Most implementations invoke the 
adiabatic approximation, by assuming that the exchange- 
correlation (XC) kernel, the second functional deriva- 
tive of the XC energy with respect to density, is fre- 
quency independent. Maitra and coworkers [65] showed 
that a frequency-dependent XC kernel is necessary to 
correctly describe a state with a strong double excita- 
tion mixing. Double excitation energies using the adi- 
abatic form of the quadratic response are simply the 
sums of two single excitation energies, 0, [H, [HJ unlike 
early claims to the contrary. [68| These trivial double ex- 
citation energies may not be found within the Tamm- 
Dancoff approximation. ^ Maitra et al. also proposed 
dressed TDDFT [65] to remedy this deficiency; deriving 
a frequency-dependent XC kernel in which single exci- 
tations are mixed with spectrally isolated double excita- 
tions. This approach has been applied and tested. [70 ilSj 
Its major weakness is the need to assign the mixed single 
and double excitations a priori. 

Many additional density functional methods have been 
proposed for doubly-excited systems. Casida had added 
non-DFT many-body polarization propagator correc- 
tion based on the Bethe-Salpeter equation to the XC 
kernel. [54] This nonadiabatic XC kernel includes the 
dressed TDDFT kernel as a special case if the ground 
state has a closed shell. Additional frequency-dependent 
XC kernels accounting for double excitations in fi- 
nite and correlated systems were derived recently from 
the Bethe-Sal2eter equation, jil Q The spin-fiip (SF) 
approach |75l - [79| is also a promising method for generat- 
ing double excited states from a triplet single reference 
state, in which electrons are excited to an orbital with 
a different spin. The SF approach can be employed to 
both wave-function- [75|-[78| and TDDFT [79| based meth- 
ods. In the original implementation of SF-TDDFT, the 
coupling of SF excitations entered the linear response 
equation only through the Hartree-Fock exchange com- 
ponent in the hybrid XC functionals. Wang and Ziegler 
overcame this limitation by a noncollinear formulation 
of the XC potential. [80i] Recent assessment and appli- 
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FIG. 1: Schematics of self-consistent calculation of 
core-excited states. The EGA (left) replaces a core with 
the next highest element in the periodic table, the FGH 
(middle) fixes the occupation of the core-orbital, 
neglecting the excited electron, whereas the XGH 
(right) fixes the occupation of both orbitals. 



cations of SF-TDDFT are given in Refs. ISlHSJ. The 
real-time approach has also been used to obtain double 
excitations. (sBl Iscj A time-independent DFT method for 
multiple excitations was proposed recently, [st", '88^ In this 
approach multiple excited states were obtained from an 
optimized effective potential (GEP) eigenvalue equation 
with orthogonality constraints. Gomputing of the GEP of 
a polyatomic molecule, however, is not easy, and there are 
still many open questions in the DFT response formal- 
ism for excited states. Time-dependent density matrix 
functional theory (TDDMFT), proposed by Giesbertz et 
al.,[89] may also account for double excitations. Adopt- 
ing these methods to better describe the processes that 
contribute to various nonlinear spectroscopy techniques, 
and design experiments which test the assumptions un- 
derlying them, is an ongoing open challenge. 



III. THEORY 

A. Survey of Approximations for Single Core Holes 

The three most common approximate descriptions of 
core excitations are illustrated in Fig. [TJ The equivalent 
core hole approximation (EGA), also known as the Z 
approximation, [90, 91] replaces the core hole by an addi- 
tional nuclear charge, and was employed in our previous 
simulation studies of X-ray nonlinear spectroscopy sig- 
nals. [92-96] The EGA can be easily implemented within 
the routine features of standard quantum chemistry pack- 
ages, and easily represents DGESs using an extra two 
nuclear charges. However, it is only applicable to deep 
core-holes as it neglects the effect of the chemical envi- 
ronment, changes the spin symmetry of a single core hole 



state [97] and cannot account for core hole delocalization 
and migration. 

The full core hole (FGH) approximation [98| improves 
on the basic EGA model by employing orbitals deter- 
mined self-consistently using the fixed core hole configu- 
ration as shown in the middle panel of Fig. [TJ The di- 
rect static exchange (STEX) model [99|, llQOj is one type 
of FGH with the occupied Hart ree-F ock orbitals and the 
improved virtual orbitals (IVO) |lQl( of the ionic {N — 1)- 
electron system. Unlike the EGA, STEX includes orbital 
relaxation and can be applied to both deep and shal- 
low core- holes. Interactions between the excited electron 
and the other — 1 inner electrons, however are still 
neglected. We have used STEX to study linear and non- 
linear X-ray spectroscopy signals [33|, llQ2[ llQ3| in small 
organic molecules. A detailed analysis of ST EX an d com- 
parison with REW-TDDFT is given in Refs. [Tq2| and [33l. 
Gurrent implementations of STEX are limited to single 
core hole states. Spin coupling of the two core electrons 
compHcates the DGES wave-function, making it hard 
to map to an effective single particle Hamiltonian like 
that in STEX. Self-consistent field (SGF) calculations of 
DGESs are also numerically tricky, and frequently fail to 
converge. Previous calculations of do uble core hole states 
using GASSGF, [13, 39, 104] MRGL [ToBI . [ToGj and ADG 
[lQ5l llQ7l4lQ9| are expensive and had only been applied 
to small molecules such as NH3 and GH4. R elatiyistic 
corrections may be necessary for core electrons. llQ5l fTTol 
The much cheaper ASGF method, [Joi [TToHml like 
STEX, suffers fr om the SGF convergence problem for 
core hole states |ll3| . Moreover, running many SGF 
calculations with combinations of core-holes and excited 
electrons is tedious. 

We employ a third approach, REW-TDDFT, first pro- 
posed by Stener et al. in 2003 [26] and further developed 
in additional studies. [27., ^] REW-TDDFT only consid- 
ers electrons excited from a defined set of relevant or- 
bitals (the restricted excitation window), allowing to ob- 
tain high-lying states in the excitation spectrum without 
calculating the lower states. A similar restricted chan- 
nel approach was suggested in Ref. without orbital 
relaxation in the field of the core hole. Like TDDFT, 
the complex polarization propagator method [114] is also 
based on response theory. With this method, the absorp- 
tion of the system at a given frequency can be calculated 
by solving a response matrix equation; the interesting 
energy region can be sampled without solving explicitly 
for the excited states. This method was applied in Ref. 

We have recently found [33| REW-TDDFT to be more 
accurate in predicting frequency splitting in XANES 
and computationally less demanding than STEX. Brena 
et al. had pointed out that TDDFT core excita- 
tion energies have larger absolute errors than those 
from STEX or transition-state calculations, because 
TDDFT does not account for orbital relaxation and self- 
interaction of core electrons. |ll6( We found the same 
trend. [33] Apart from an overall shift, TDDFT core ex- 
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citation energies agree w ith the XANES of many sys- 
tems very wen.[3l],[32|,[Ii7i] The TDDFT core-edge energy 
can be improved by apply ing_the Per dew- Zunger self- 
interaction correction scheme, [l 18', "ll9] or by employing 
rece ntly-developed core- valence hybrid functional, [120'- 
Il22| long-range corrected hybrid functional with short- 
range Gaussian attenuation [123] and short-range cor- 
rected functionals.[29, 30] These techniques yield core- 
edge energies with errors less than leV on a series of 
small molecules. (291 [301] More precise benchmarks will be 
needed to compare the accuracy of STEX and REW- 
TDDFT. 




FIG. 2: The XDQC technique. 



B. Double-core Excitations 

Adiabatic TDDFT cannot deal with double 
excitations. [65| Practical frequency-dependent XC 
kernels for medium or large molecular systems have yet 
to be developed. Both TDDFT and the excited state 
core hole (XCH) method [124] treat singly core-excited 
state (SCES) very well. The excited electron in XCH 
approximation is included self-consistently through a 
full core hole state. XCH has been applied to solid 
state X-ray absorption spectra with pseudopotentials, 
and applying it to molecular system poses no additional 
difficulty. We will use XCH in this paper. Pseudopo- 
tentials can be constructed for single [125] and double 
core hole states. An effective potential which ignores 
the polarizability of the core hole in different chemical 
environments will introduce errors into practical calcu- 
lations. A fixed pseudopential of a core hole, like the 
EGA, ignores all core hole dynamics. 

Here we employ an approach that combines XCH and 
REW-TDDFT. We ffist run an SCF calculation to get a 
reference state SCESO, with a core hole and an excited 
electron. This reference is then used to run a REW- 
TDDFT calculation to obtain excited states with two 
core-holes and two excited electrons. Adiabatic TDDFT 
only treats single excitations of the SCESO reference, not 
all double core excited states are accounted for. Addi- 
tional DCESs can be found by starting REW-TDDFT 
at different SCES references, obtained by permuting the 
occupied and virt ual o rbitals of SCESO (orbital approxi- 
mation as in Ref. Il02l ). 

Unrestricted reference-based TDD FT is known to suf- 
fer from spin contamination. Il26( A truncated rank of 
excitations from one component of a spin multiplet gener- 
ates an incomplete configuration space for the total spin 
operator 6'^. [57i , il27. . J28] Stated differently, higher rank 
excited configurations are needed to represent a pure 
spin state. Our current TDDFT calculations on a spin 
symmetry-broken reference state are no exception, but 
since the singlet-triplet energy splitting is negligible com- 
pared to the large energies of core excitations, spin con- 
tamination does not strongly affect the calculated core 
excitation energies. Casida and coworkers suggested to 
use the difference of the total spin between excited states 



and the reference state (A ( 5^)) in order to ffiter out un- 
physical excited states. jl26( We follow his scheme by only 
including excited states with small A(S'^)s in our calcu- 
lations. 

The two sequential core excitations cannot be treated 
in the same fashion. After the initial excitation, the va- 
lence electrons relax self-consistently (XCH) in the field 
of the core hole, and the second excitation is treated using 
response theory (REW-TDDFT). This treatment uses 
two potentials to perturb the valence band during the 
second excitation: a strong, long duration core hole po- 
tential which must be treated nonperturbatively, and the 
ultrafast second photon excitation, which can be trun- 
cated at ffist order using response theory. This physi- 
cally intuitive order in which relaxation precedes photon 
excitation strongly effects the signal, as we show later 
by reversing it, treating the original ffist core excitation 
using REW-TDDFT and the original second using XCH. 



IV. THE DQC SIGNAL 

The DQC signal employs four pulses with wave vectors 
satisfying kni = ki+k2— ks, and the transmission change 
of the kin pulse is recorded versus the three interpulse 
delays (t 1,^2, ^3) (as shown in Fig. [2]). The applied electric 
field is 
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Ej{t-Tj)^C.C.. 



(1) 



We represent the pulses as 

t-rj)= €j{t - fj) exp ' r - iujj(t - fj)] (2) 

with central frequencies cjj, wave vectors kj (k4 = km), 
and Gaussian envelopes 



ej(t-Tj) 



E^e^p [-(t-f,)V2a|] 



(3) 



with amplitudes E^^ temporal widths aj and envelope 
centers fj. The XDQC signal depends on the three 
(positive) delay times {tj = fj+i — fj) or their Fourier 
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conjugates (l^i, fl2^ ^s)- Invoking the rotating wave 
approximation and assuming temporally well-separated 

I 



pulses, [l29^ th e sig nal is given by two the loop diagrams 
of Fig. [3] tllO, 



Siu{Qs, 1^2, ^l) = ^III;a(^3, ^2, ^l) + ^III;b(^3, ^2, ^l), 



(4) 



where 



^ (1^3 - CJe'g + he'g){^2 - f g + ^7/^)(^l - ^eg + i^feg) 



and 



^lll-B^^^ (^3, ^2, ^l) 



E 

fe'e 



e'^))(/^/'e^2(^2 



(6) 



r 



Here, is the spectral envelope of the j pulse (given 

by the Fourier transform of Eq. [3j), vi,,a are the tensor 
components of the transition dipole /i^s, ^rs is the tran- 
sition frequency between the states r and s, and is 
phenomenological parameter describing the inverse life- 
time of the core excited state. Different experimental 
techniques measure various projections of the full three- 
dimensional response in Eq. IH Since XDQC resonances 
show up along 0.2 we shall display two 2D projections of 
this 3D signal 

/oo 
e-i03*3 5jjj(Q3^Q2,Oi)dn3- (7) 
-00 

where we hold fixed, and 

/oo 
e-*"i*i5iii(1^3,f^2,f^i)dl^i (8) 
-00 

for a fixed ti. Fixed, nonzero ^1^3 avoid the possibility 
of overlapping pulses contributing to the signal. There 
are three independent tensor components of the signal in 
isotropic media which depend on contractions over differ- 
ent field polarization vectors. The signal in an isotropic 
sample is a linear combination of various contractions 
over the tensor components vi,,a in Eq. |4j The rotation- 
ally averaged signal with all parallel pulse polarizations 
is given by Eq. | 



V. COMPUTATIONAL DETAILS 

The optimized geometry of formamide was taken from 
Ref. m. XCH calculations were performed by converg- 
ing the electron configuration with a designated core hole 
and an excited electron on the lowest unoccupied molec- 
ular orbital (LUMO). REW-TDDFT calculations were 



performed with the Tamm-Dancoff approximation. The 
calculation of DCESs was described in Sec. IIIIBI First a 
set of XCH-relaxed orbitals are acquired, and the REW- 
TDDFT equations are solved in this basis. Transition 
dipole matrix elements between singly and doubly ex- 
cited states are evaluated between Kohn-Sham determi- 
nants with nonorthogonal orbitals using Eq. 10 from Ref. 
33: 



config. 



(9) 

where ^a,b are the SCES and DCES wavefunction re- 
spectively, d is the transition dipole operator, am and hn 
are configuration interaction (CI) coefficients for differ- 
ent excited configurations (m and n) of the SCES A and 
DCES B, respectively. 



(10) 



is the transition dipole matrix between single excitation 
configurations m and n, Cip^rn.A and Cjq^n,^ are MO coef- 
ficients for the configurations m and n of the SCES and 
DCES, respectively. 



nmn 



k,l 



(11) 



is the overlap matrix between the MOs of the config- 
urations m and n of state A and B, in Eqs. [TOl 
and [TTl are basis functions and i^j^p^q^k^l are indices for 
these basis functions. Minor (S^"^)ij denotes the (z,j) 
minor of the matrix S^^. Al l calculations were car- 
ried out at the B3LYP[132l [i3l/cc-pVTZ[l^ level us- 
ing a modified vers ion of the quantum chemistry package 
NWChem.||i3,[i3i 
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FIG. 3: (left) The two diagrams contributing to the double quantum coherence signal, (right) Molecular structure 
and REW-TDDFT level scheme. \OjNi) refers to the double core-excited states with the 01s electron excited to the 
jth virtual orbital of the self-consistent core state, while the Nls to ith orbital excitation is obtained through 

response theory, as described in the text. 
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FIG. 4: Calculated (grey) nitrogen (left) and oxygen 
(right) K-edg e XA NES for formamide. Experimental 
EELS spectra |l4Q( are given as black lines, and power 

spectra of the pulses used in the calculation of the 
2D-QCS signals as dashed lines. The simulated energies 
were shifted (+13.065 eV for nitrogen and +14.5 eV for 
oxygen K-edge) to fit the EELS signals [ 140] 



VI. SIMULATION RESULTS 
A. XANES 

The SCES frequencies (eg) and transition dipoles (fieg) 
were obtained using the REW-TDDFT response ap- 
proach and a ground state reference wavefunction as de- 
scribed in Sec. [Ill Te, the lifetime broadening of the core- 
excited state |e), is set to 0.1 eV. The XANES signal 

5'XANES(C^) = - 7 ^""^^^2 ^, 2 • (^^) 

is shown in Fig. |4j X-ray and electron scattering in the 
core energy range are known to resemble each other in 
the gas phase. [136] The experimental electron energy loss 
spectra (EELS) are shown for comparison. 



The two lowest-frequency peaks in the experimental 
nitrogen K-edge EELS are split by 0.95 eV, and the low 
energy red peak has an intensity x 0.708 relative to that 
of the blue peak. Our REW-TDDFT XANES simula- 
tions show a larger, 1.29eV splitting, and a more intense 
blue component (x2.82 of the red peak). The simulated 
oxygen edge spectrum is much closer to experiment, and 
the splitting between the lowest energy peak and higher- 
energy transitions is reproduced. We see a strong, intense 
core-edge distinct from the higher energy shoulder. Ex- 
perimental peaks significantly higher than the ionization 
energy are broadened by coupling to unbound photoelec- 
tron states, a process not included in our simulations. 



B. The XDQC Signal 

The Siu{ts = 5fs,(^2,^i) signal is displayed in Figs. 
[5] and [H In this plot, resonances on the Qi axis reveal 
SCES {oOeg)^ and along the Q2 axis we see the DCES 
(ojfg) frequencies. Time-evolution during the ts period 
results in a phase depending on the final state |e'); these 
phases are displayed using the colors for each peak in the 
2D stick spectrum in the top row. 

In Fig. [5l we excite the oxygen core electron first, cre- 
ating resonances along l^i, and monitor their infiuence on 
the nitrogen core excitations, in the DCESs along the Q2 
axis. This pulse order is marked ONNO. In the NONO 
spectra shown in Fig. [6] where nitrogen core excitation is 
first, we see the infiuence of nitrogen core excitations on 
the oxygen core excitations. The ONNO signal has reso- 
nances inQi and Qs representing core-excitation with the 
same element, oxygen, in both diagram A {Q^ ^ ^e'g)i 
and diagram B (Q.^ ^ Ufe^). Fig. [5] shows a typical pat- 
tern of the XDQC signal, a series of features lying on 
parallel lines of roughly constant 1^2 — ^1 (diagonal char- 
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acter). This pattern indicates weak correlation between 
the DCESs and SCESs: the energies of the DCESs are 
roughly the sum of the energies of two SCESs. The oxy- 
gen core-excitations do not substantially affect the nitro- 
gen core-excitation spectra. In contrast, the NONO sig- 
nal (Fig. [6]) shows a more complicated variation with l^i, 
in which the diagonal pattern is blurred. Two |e) reso- 
nances lie below the diagonal at 407 eV, suggesting 
that self-consistent relaxation in the field of the nitro- 
gen core hole after the first pulse leads to a set of DCESs 
during t2 in which the oxygen and nitrogen excited states 
strongly interfere. To more closely examine the difference 
between these two pulse orders, we compare the peak 
splittings between the XANES and 2D spectra. The en- 
ergy difference between A' and A is around 1.29 eV in 
XANES. A' and A correspond to the two strong features 
in Fig. [5] (shown more clearly in the left panel of Fig. [7]), 
whose energy difference is 1.34 eV, indicating that the 
oxygen core hole does not appreciably affect the accessi- 
ble virtual orbitals during nitrogen single core excitation. 
However, the energy difference of the two strong features 
on the Q2 axis in Fig. [Gfshown more clearly in the right 
panel of Fig. [7j), shows that the dominant strong oxygen 
single core excitation is shifted differently by the vari- 
ous nitrogen single core excitations. Comparison of the 
XDQC to the XANES spectra provides useful insights on 
correlations between specific core excitations. 

To better illustrate the relationship between the 2D- 
XDQC signal and the linear absorption, we display the 
ONNO and NONO signals with the XANES marginals 
in Fig. El We shift the Nls and 01s XANES on the 
Q2 axis to match the strongest features. If the core ex- 
citations were uncoupled, the peaks in for a fixed 
fti should reproduce the XANES of the core resonant 
with the second pulse, with one peak quenched due to 
Pauli repulsion with the previously excited core-electron. 
This might explain the strong diagonal character of the 
ONNO signal, in which the strong peak A is repeated 
for increasing Qi. This interpretation does not hold for 
the NONO signal, with its group of four strong peaks 
around l^i = 405eV and ft2 = 950eV. The pattern does 
not match a shifted XANES signal, implying strong cor- 
relations in the DCESs. Comparing the linear absorption 
(Fig. SI) to the ONNO XDQC signal for ts =5fs (Fig. [5]), 
we find one medium and two stronger features in the ni- 
trogen K-edge XANES, and one strong, two weak and 
one medium intense features for oxygen K-edge XANES. 
The corresponding peaks are marked as A', A, B , C, E, 
F, D in Fig. H The N and O K-edge XANES are well 
described by a single particle picture. The dominant par- 
ticle MOs in the CI expansions of these SCESs are shown 
on Fig. E 

In the ideal zero-DCES-SCES correlation case we ex- 
pect to see two parallel Hues of strong features in Fig. [5] 
corresponding to A and B in the nitrogen XANES, but 
we only see one clear set of "diagonal" resonances. The 
peaks corresponding to states E and F in the oxygen 
XANES are much weaker than that of D, but their cor- 



Total 



Diagram A 



Diagram B 




525 530 535 540 545 550 525 530 535 540 545 550 525 530 535 540 545 550 

njey QJey QJey 

FIG. 5: The Sm{t3 = 5fs,172,l^i) ONNO signal with 
XXXX polarization configuration. The total signal (left 
column) is the sum of the contributions from diagram A 
(middle column) and diagram B (right column) of Fig. 

[3l Each circle in the stick spectra (top row) has a 
complex contribution to the signal from a combination 
of states, with the radius of the circle proportional to 
the square root of the amplitude, and colored according 

to the phase of the contributing peak. The following 
three rows show the absolute value, real and imaginary 
parts of the complex signal after convoluting with a 
Lorentzian of width O.leV. All signals were scaled so 
that abs(6'jfj*(t3 = 5fs, 1^2,^1)) has a maximum value of 

one. 



responding peaks in the XDQC spectra (in Fig. [5]) are as 
strong as as that of D. State E shares the same particle 
MO as that of A\ implying a strong coupling between E 
and A', which explains the corresponding strong XDQC 
signal. But for the XDQC peak corresponding to F, it is 
very difficult to compare the coupling strength of F and 
A' to that of D and A\ simply by their dominant particle 
MOs. Even the calculated energy of the DCES is close 
to the sum of two single excitations, the wavefunction 
is not the simple outer product of two SCES wavefunc- 
tions. The single particle picture can not explain why 
the XDQC peak corresponding to F also becomes strong, 
suggesting that many-body effects dominate the signal. 

As discussed in Sec. [Ill the DCES (|/)) strongly de- 
pend on the order of core excitation by the first two 
pulses, showing that XDQC is sensitive to the order of 
the formation of doubly core-excited states. Our XDQC 
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Total 



Diagram A 



Diagram B 




395 400 405 410 415 420 395 400 405 410 415 420 395 400 405 410 415 420 

njey Qi/ev njey 

FIG. 6: same as Fig. O but for the pulse sequence 
NONO. 



simulations also show large differences between the two 
protocols explained in the following, reflecting the ap- 
proximate treatment of many-body effects. In Figs. [5] 
and [6] we assumed the |/) states were relaxed in the field 
of the first (core hole a), and calculated using response 
theory with respect to the second excited core (core hole 
b, protocol i). In Fig. [9] we present an example of the op- 
posite protocol for the ONNO technique: the |/) states 
are generated using SCF relaxation for core hole b and 
then response theory for core hole a (protocol ii). We 
observe a different blurred line pattern in the spectra ob- 
tained with protocol ii (see Fig. [9|). Different levels of 
theory show a different type of DCES-SCES correlation. 
The protocol i diagram A signal also shows a stronger 
diagonal component than the corresponding protocol ii 
signal, which displays a different density of two-particle 
states on the 1^2 axis. Comparing XDQC signals from dif- 
ferent pulse orders with experiments should reveal how 
well those theories treat the very specific correlation be- 
tween a DCES and SCES. This could be valuable in fu- 
ture methodology development for DCESs. 

Another portion of the signal, 5'iii(1^3, 1^2, = 5fs), is 
shown in Figs. [10] and [TTJ In these graphs, resonances 
along Q2 correspond to doubly excited states (ojfg)^ and 
fts shows the final single excited states (c^e'^ and cj/eO 
in the XDQC process. As with the Siu{ts = 5fs, 1^2,^1) 
signals, there are diagonal characteristics in the ONNO 
spectra (Fig. [TQ|) which are absent in the NONO spectra 
(Fig. [TT]) . This is further evidence for the increased cor- 



relations for the states probed by the ONNO process. In 
the 5'iii(t3 = 5fs,r^2,^i) signal the resonances for both 
A and B diagrams are the same, only the phases for two 
contributions differ. The 5'iii(l]3, 1^2, ^1 = 5fs) signal has 
qualitatively different peaks along for the two dia- 
grams, as expected, since these diagrams differ in this 
time period. 



VII. CONCLUSIONS 

Typically the overlaps of valence holes or electron or- 
bit als are higher than those of core holes or electrons. 
The correlations between valence excitations are stronger 
than those between core excitations, which makes the 
corresponding DQC signals more complicated than those 
in the present work.[137j We have demonstrated qual- 
itative differences in the stimulated XDQC signals de- 
pending on the order of core-excitations at different ele- 
ments in a small organic molecule. As an experimental 
technique, XDQC may be used to fingerprint different 
theoretical approaches for modeling core-excitation, by 
dissecting dynamical doubly-core-excited resonances cor- 
responding to the effects of valence relaxation induced by 
the core hole. Additional work will be required to fur- 
ther analyze the spectra reported here. First, a detailed 
description of the valence dynamics, and a comparison 
of the TDDFT theory with the more easily interpreted 
orbital theory should highlight the effect of electron cor- 
relation on the XDQC signal. Second, with a higher level 
of theory to describe DCESs, the variation of the Q2 res- 
onances with ti can be used to measure how the order of 
core-excitation affects the two-particle density of states. 
The core hole has a unit charge which creates a very 
strong field over atomic length scales. If the delay be- 
tween the first two pulses is varied, the changes in the 
I/) resonances along Q2 will reflect the order of pertur- 
bation theory necessary to correctly model this strong 
interaction with the valence band. For short ti, linear 
response theory alone, or supplemented with a full or 
partial core hole orbital transformation may adequately 
capture this effect, as it does in linear absorption. We 
demonstrated significant differences between the XDQC 
signals predicted by the calculation protocols. Interpola- 
tion between the extremes of linear response (first or- 
der), and a full SCF calculation will require a higher 
level of theory to describe many-body effects in the core- 
excitation. 

In the simulations presented here the relaxation of sin- 
gle and double core-excited states, is treated phenomeno- 
logically. We only include population decay through 
Auger or a radiative process that destroys the core hole. 
We further use the same decay rate for all SCESs and 
DCESs, independent of the orbitals involved. Interac- 
tion of the system with a bath, composed of vibrational 
and valence-electron degrees of freedom, could introduce 
additional pure dephasing. Nonlinear spectroscopy has 
been successful in lower frequency regimes (NMR to the 
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I I I III / II III/ 

J \ I I \ \ I I I U—i I \ I I \ \ I 1 I I I I ^ 

525 530 535 540 545 550 395 400 405 410 415 420 



FIG. 7: Comparison of the absolute parts of the ONNO (left) and NONO (right) Siu{ts = 5fs, 1^2,^1) signals. 

XANES spectra are shown in the marginals. 




FIG. 8: Dominant MOs of single particle orbitals of different SCESs discussed in Section IVl B[ (a) For peak E and 

A^ (b) For peak F (c) For peak D. Peaks are labeled in Fig. |4j 



visible) at disentangling different und erlyin g mechanisms 
for excited state decay and dephasing.[l38] Similar effects 
are expected in core-excitation spectra. Estimating their 
magnitude, comparing them to population decay rates, 
and designing experiments that distinguish between them 
is an interesting future topic. One manifestation of these 
effects in existing X-ray experiments are fluorescence vs. 
Raman signals, which are controlled by the ratio of pure 
dephasing and population relaxation rates. 
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Appendix A: Orientational Tensor Averaging 



The dipole and electric field polarization tensor contri- 
butions to the signal must be rotationally averaged over 
all possible orientations of the molecule. The rotationally 
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Diagram A 



Diagram B 



Total 



Diagram A 



Diagram B 



> 955 



C 950 



> 955 
V CU 



G 950 



i i><^ 






' A 

1 - 




. 1 








111 « 1 * 







525 530 535 540 545 550 525 530 535 540 545 550 

FIG. 9: Comparison of the two protocols to calculate 
the double-core excited states |/) for the all-parallel 
ONNO Siu{ts = 5fs, 1^2,^1) signal as discussed in Sec. 
IVIBI Contributions from diagram A (left column) and 
B (right column) with protocols i (top row) or protocol 
ii (bottom row). All graphs were multiplied by the same 
scaling factor used in Fig. [5] 



Total 



Diagram A 



Diagram B 




940 945 950 955 960 965 940 945 950 955 960 965 940 945 950 955 960 965 

^s/eV /eV njey 

FIG. 10: 8111(^3^2, h = 5fs) for the ONNO pulse 
configuration. 




940 945 950 955 960 965 940 945 950 955 960 965 940 945 950 955 960 965 

/eV n^/ey 

FIG. 11: (1^3, ^2,^1 = 5fs) for the NONO pulse 
configuration. 



averaged 1 139j signal is 



crot ^ 

^III — ^ 




4 


-1 


-1 


1 


4 


-1 


1 


-1 


4 



cos 612 COS 6>34 
COS Ols COS O24 
COS 6>i4 COS O23 

(Al) 

where cos Oij = • ej , is the the polarization of the ith 
pulse polarization vector, and Sa,/3,j represent various 
contractions over the tensor components of the response; 



^III — / . ^] 



V'2 V2 ^1 ^1 
III ' 



On 



1^2 ^1 ^1 ^\ 
III ' 



and 



'^iii — ^iii 



(A2) 



(A3) 



(A4) 



For the all parallel pulse configuration, the signal is pro- 
portional to the sum of all three 



11 



[1] T. Popmintchev, M.-C. Chen, P. Arpin, 
M. M. Murnane, and H. C. Kapteyn, 

Nat. Photonics 4, 822 (2010) 

[2] F. Krausz and M. Ivanov, Rev. Mod. Phys. 81, 163 
(2009). 

[3] P. Emma, R. Akre, J. Arthur, R. Bionta, C. Bost- 
edt, J. Bozek, A. Brachmann, P. Bucksbaum, R. Cof- 
fee, F.-J. Decker, Y. Ding, D. Dowell, S. Edstrom, 

A. Fisher, J. Frisch, S. Gilevich, J. Hastings, G. Hays, 
P. Hering, Z. Huang, R. Iverson, H. Loos, M. Messer- 
schmidt, A. Miahnahri, S. Moeher, H.-D. Nuhn, G. Pile, 

D. Ratner, J. Rzepiela, D. Schultz, T. Smith, P. Stefan, 
H. Tompkins, J. Turner, J. Welch, W. White, J. Wu, 
G. Yocky, and J. Galayda, Nat. Photonics 4, 641 (2010). 

[4] J. Ullrich, A. Rudenko, and R. Moshammer, Ann. Rev. 
Phys. Chem. 63, 635 (2012). 

[5] L. Young, E. P. Kanter, B. Krassig, Y. Li, A. M. 
March, S. T. Pratt, R. Santra, S. H. Southworth, 
N. Rohringer, L. F. DiMauro, G. Doumy, C. A. Roedig, 
N. Berrah, L. Fang, M. Hoener, P. H. Bucksbaum, 
J. P. Cryan, S. Ghimire, J. M. Glownia, D. A. Reis, 
J. D. Bozek, C. Bostedt, and M. Messerschmidt, 
Nature 466, 56 (2010), 

[6] M. Hoener, L. Fang, O. Kornilov, O. Gessner, S. T. 
Pratt, M. Giihr, E. P. Kanter, C. Blaga, C. Bostedt, 
J. D. Bozek, P. H. Bucksbaum, C. Buth, M. Chen, 
R. Coffee, J. Cryan, L. DiMauro, M. Glownia, E. Hosier, 

E. Kukk, S. R. Leone, B. McFarland, M. Messerschmidt, 

B. Murphy, V. Petrovic, D. Rolles, and N. Berrah, 
Phys. Rev. Lett. 104, 253002 (2010). 

[7] P. Elliott, S. Goldson, C. Canahui, and N. T. Maitra, 

Chem. Phys. 391, 110 (2011). 
[8] J.-P. Connerade, Highly Excited Atoms (Cambridge 

University Press, 1998). 
[9] R. Santra, N. V. Kryzhevoi, and L. S. Cederbaum, 

Phys. Rev. Lett. 103, 013002 (2009). 
[10] K. Ueda and O. Takahashi, J. Electron. Spectrosc. Re- 

lat. Phenom. 185, 301 (2012). 
[11] N. V. Kryzhevoi, R. Santra, and L. S. Cederbaum, 

|J. Chem. Phys. 135, 084302 (2011), 
[12] P. Salen, P. van der Meulen, H. T. Schmidt, R. D. 

Thomas, M. Larsson, R. Feifel, M. N. Piancastelli, 

L. Fang, B. Murphy, T. Osipov, N. Berrah, E. Kukk, 

K. Ueda, J. D. Bozek, C. Bostedt, S. Wada, R. Richter, 

V. Feyer, and K. C. Prince, Phys. Rev. Lett. 108, 

153003 (2012). 

[13] M. Tashiro, K. Ueda, and M. Ehara, Chem. Phys. Lett. 

521, 45 (2012). 
[14] J. Kim, S. Mukamel, and G. D. Scholes, Acc. Chem. 

Res. 42, 1375 (2009). 
[15] S. Mukamel, R. Oszwaldowski, and D. Abramavicius, 

Phys. Rev. B 75, 245305 (2007). 
[16] W. Zhuang, D. Abramavicius, and S. Mukamel, Proc. 

Natl. Acad. Sci. 102, 7443 (2005). 
[17] A. Nemeth, F. Milota, T. Mancal, T. Pullerits, J. Sper- 
ling, J. Hauer, H. F. Kauffmann, and N. Christensson, 

J. Chem. Phys. 133, 094505 (2010). 
[18] S. Mukamel, R. Oszwaldowski, and L. Yang, 

[jTChem. Phys. 127, 221105 (2007) 
[19] D. Karaiskaj, A. D. Bristow, L. Yang, X. Dai, 

R. P. Mirin, S. Mukamel, and S. T. Cundiff, 



[2o; 

[21 



[22" 
[23^ 

[24 



[25; 
[26; 
[27; 
[28; 

[29; 
[3o; 

[31 

[32; 

[33 

[34; 
[35; 
[36; 

[37" 

[38; 

[39; 

[4o; 

[41 

[42; 

[43; 
[44; 
[45; 

[46; 
[47; 



|Phys. Rev. Lett. 10 4, 117401 (2010)' 



K. W. Stone, K. Gundogdu, D. B. Turner, X. Li, S. T . 
Cundiff, and K. A. Nelson, S cience 324, 1 169 (2009)^ 
X. Dai, M. Richter, H. Li, A. D. Bristow, 
C. Falvo, S. Mukamel, and S. T. Cundiff, 



Phys. Rev. Lett. 108, 193201 (2012) 

J. Orenstein, Phys. Today 65, 44 (2012). 

G. D. Mahan, Condensed Matter in a Nutshell (Prince- 
ton University Press, 2011). 

D. Fausti, R. L Tobey, N. Dean, S. Kaiser, A. Dienst, 
M. C. Hoffmann, S. Pyon, T. Takayama, H. Takagi, and 
A. Cavalleri, Science 331, 189 (2011). 
L V. Schweigert and S. Mukamel, Phys. Rev. A 78, 
052509 (2008). 

M. Stener, G. Fronzoni, and M. de Simone, Chem. 
Phys. Lett. 373, 115 (2003). 

N. A. Besley and A. Noble, J. Phys. Chem. Gill, 3333 
(2007). 

S. DeBeer-George, T. Petrenko, and F. Neese, J. Phys. 
Chem. A 112, 12936 (2008). 

N. A. Besley, M. J. G. Peach, and D. J. Tozer, Phys. 
Chem. Chem. Phys. 11, 10350 (2009). 
N. A. Besley and F. A. Asmuruf, Phys. Chem. Chem. 
Phys. 12, 12024 (2010). 

W. Liang, S. A. Fischer, M. J. Frisch, and X. Li, J. 

Chem. Theory Comput. 7, 3540 (2011). 

K. Lopata, B. E. V. Kuiken, M. Khalil, and N. Govind, 

J. Chem. Theory Comput. 8, 3284 (2012). 

Y. Zhang, J. D. Biggs, D. Healion, N. Govind, and 

S. Mukamel, J. Chem. Phys. 137, 194306 (2012). 

N. Ottosson, E. F. Aziz, H. Bergersen, W. Poka- 

panich, G. Ohrwall, S. Svensson, W. Eberhardt, and 

O. Bjorneholm, J. Phys. Chem. B 112, 16642 (2008). 

N. Ottosson, E. F. Aziz, L L. Bradeanu, S. Legendre, 

G. Ohrwall, S. Svensson, O. Bjorneholm, and W. Eber- 
hardt, Chem. Phys. Lett. 460, 540 (2008). 

H. Ikeura-Sekiguchi, T. Sekiguchi, Y. Kitajima, and 

Y. Baba, App. Surf. Sci. 169-170, 282 (2001). 

D. P. Chong, J. Electron. Spectrosc. Relat. Phenom. 184, 16 4 (2011 

H. Agren and H. J. A. Jensen, Chem. Phys. 172, 45 

(1993). 

O. Takahashi, M. Tashiro, M. Ehara, K. Yamasaki, and 
K. Ueda, Chem. Phys. 384, 28 (2011). 

jL Jadzyn and Swiergiel, 

iPhys. Chem. Chem. Phys. 14, 3170 (2012)1 
Y. L. Li, L. Huang, R. J. Dwayne Miller, T. Hasegawa, 
and Y. Tanimura, J. Chem. Phys. 128 , 234507 (2008) 
A. Paarmann, M. Lima, R. Chelh, V. V. Volkov, 
R. Righini, and R. J. D. Miller, Phys. Chem. Chem. 
Phys. 13, 11351 (2011). 

J. Jortner and S. A. Rice, J. Chem. Phys. 44, 3364 
(1966). 

D. Sangalli, P. Romaniello, G. Onida, and A. Marini, 
J. Chem. Phys. 134, 034115 (2011). 
L. Serrano- Andres, M. Merchan, L Nebot-Gil, R. Lindh, 
and B. O. Roos, J. Chem. Phys. 98, 3151 (1993). 
K.N. amd Haruyuki Nakano and K. Hirao, Int. J. Quan- 
tum Chem. 66, 157 (1998). 

C. Hsu, S. Hirata, and M. Head- Gordon, J. Phys. 
Chem. A 105, 451 (2001). 



12 



[48] J. H. Starcke, M. Wormit, J. Schirmer, and A. Dreuw, [77 

Chem. Phys. 329, 39 (2006). 
[49] J. H. Starcke, M. Wormit, and A. Dreuw, J. Chem. [78 

Phys. 131, 144311 (2009). 
[50] S. Knippenberg, J. H. Starcke, M. Wormit, and [79 

A. Dreuw, Mol. Phys. 108, 2801 (2010). 
[51] C. Angeh and M. Pastore, J. Chem. Phys. 134, 184302 [80 

(2011) . 

[52] M. Schmidt and P. Tavan, J. Chem. Phys. 136, 124309 [81 

(2012) . 

[53] E. Otero, N. Kosugi, and S. G. Urguhart, J. Chem. [82 
Phys. 131, 114313 (2009). 

[54] M. E. Casida, J. Chem. Phys. 122, 054111 (2005). [83; 

[55] M. Casida, A. Ipatov, and F. Cordova, in Time- [84 
Dependent Density Functional Theory, Lecture Notes 
in Physics, Vol. 706, edited by M. A. Marques, C. A. [85 
Ullrich, F. Nogueira, A. Rubio, K. Burke, and 
E. K. U. Gross (Spinger-Verlag Berhn Heidelberg, 2006) [86 
Chap. 16, p. 243. 

[56] P. Th. van Duijnen, S. N. Greene, and N. G. J. [87 
Richards, J. Chem. Phys. 127, 045105 (2007). 
Z. Li and W. Liu, J. Chem. Phys. 133, 064106 (2010). [ 
Y. Kurashige, H. Nakano, Y. Nakao, and K. Hirao, 
Chem. Phys. Lett. 400, 425 (2004). [89; 
V. Carravetta and H. Agren, in Computational Strate- 
gies for Spectroscopy, edited by V. Barone (John Wiley [90 
& Sons, Inc.,, 2011) p. 137205. 

O. Lehtonen, D. Sundholm, R. Send, and M. P. Jo- [91 
hansson, J. Chem. Phys. 131, 024301 (2009). 
C. M. Marian and N. Gilka, J. Chem. Theory Comput. [92 

4, 1501 (2008). 

L A. Mikhailov, S. Tafur, and A. E. Masunov, Phys. [93 
Rev. A 77, 012510 (2008). 

M. A. Marques, C. A. Ullrich, F. Nogueira, A. Rubio, [94; 
K. Burke, and E. K. U. Gross, eds.. Time- dependent 
Density Functional Theory (Springer, 2006). [95 
C. A. Ullrich, Time- dependent Density Functional The- 
ory: Concepts and Applications (Oxford University [96 
Press, USA, 2012). 

N. T. Maitra, F. Zhang, R. J. Cave, and K. Burke, J. [97 
Chem. Phys. 120, 5932 (2004). 

5. Tretiak, V. Chernyak, and S. Mukamel, Int. J. Quan- [ 
tum Chem. 70, 711 (1998). 

S. Tretiak and V. Chernyak, J. Chem. Phys. 119, 8809 
(2003). [99; 
E. K. U. Gross, J. F. Dobson, and M. Petersilka, in 
Density- Functional Theory II, Topics in Current Chem- [100 
istry, Vol. 181, edited by R. F. Nalewajski (Spinger, 
Berlin, 1996) p. 81. [101 
S. Hirata and M. Head-Gordon, Chem. Phys. Lett. 314, 
291 (1999). [102 
R. J. Cave, F. Zhang, N. T. Maitra, and K. Burke, 
Chem. Phys. Lett. 389, 39 (2004). [103; 
G. Mazur and R. Wlodarczyk, J. Comput. Chem. 30, 
811 (2009). [104; 
G. Mazur, M. Makowski, R. Wlodarczyk, and Y. Aoki, 
Int. J. Quant. Chem. Ill, 819 (2011). 

M. Huix-Rotllant, A. Ipatov, A. Rubio, and M. E. [105 
Casida, Chem. Phys. 391, 120 (2011). 

P. Romaniello, D. Sangalh, J. A. Berger, F. Sottile, [106; 
L. G. Molinari, L. Reining, and G. Onida, J. Chem. 
Phys. 130, 044108 (2009). [107 
A. I. Krylov, Chem. Phys. Lett. 338, 375 (2001). 
A. I. Krylov, Chem. Phys. Lett. 350, 522 (2001). 



Chem. Phys. 121, 12191 
and H. Agren, J. Chem. 



A. I. Krylov and C. D. Sherrill, J. Chem. Phys. 116, 
3194 (2002). 

L. V. Slipchenko and A. I. Krylov, J. Chem. Phys. 117, 
3694 (2002). 

Y. Shao, M. Head-Gordon, and A. I. Krylov, J. Chem. 
Phys. 118, 4807 (2003). 
F. Wang and T. Ziegler, J. 

(2004) . 

Z. Rinkevicius, O. Vahtras, 
Phys. 133, 114104 (2010). 

N. Minezawa and M. S. Gordon, J. Phys. Chem. A 115, 
7901 (2011). 

Z. Li and W. Liu, J. Chem. Phys. 136, 024107 (2012). 
Y. A. Bernard, Y. Shao, and A. I. Krylov, J. Chem. 
Phys. 136, 204103 (2012). 

C. M. Isborn and X. Li, J. Chem. Phys. 129, 204107 
(2008). 

W. Liang, C. M. Isborn, and X. Li, J. Chem. Phys. 
131, 204101 (2009). 

V. N. Glushkov and X. Assfeld, J. Chem. Phys. 132, 
204106 (2010). 

V. N. Glushkov and X. Assfeld, J. Comput. Chem. 33, 
2058 (2012). 

K. J. H. Giesbertz, E. J. Baerends, and O. V. Gritsenko, 
Phys. Rev. Lett. 101, 033004 (2008). 
P. Nozieres and C. T. D. Dominicis, Phys. Rev. 178, 
1097 (1969). 

W. H. E. Schwarz and R. J. Buenker, Chem. Phys. 13, 
153 (1976). 

I. V. Schweigert and S. Mukamel, Phys. Rev. Lett. 99, 
163001 (2007). 

I. V. Schweigert and S. Mukamel, Phys. Rev. A 76, 
012504 (2007). 

I. V. Schweigert and S. Mukamel, J. Chem. Phys. 128, 
184307 (2008). 

I. V. Schweigert and S. Mukamel, Phys. Rev. A 78, 
052509 (2008). 

D. M. Healion, I. V. Schweigert, and S. Mukamel, J. 
Phys. Chem. A 112, 11449 (2008). 

M. N. R. Wohlfarth and L. S. Cederbaum, J. Chem. 
Phys. 116, 8723 (2002). 

M. Cavalleri, M. Odelius, D. Nordlund, A. Nilsson, and 
L. G. Pettersson, Phys. Chem. Chem. Phys. 15, 2854 

(2005) . 

H. Agren, V. Carravetta, O. Vahtras, and L. G. M. 

Pettersson, Chem. Phys. Lett. 222, 75 (1994). 

H. Agren, V. Carravetta, O. Vahtras, and L. G. M. 

Pettersson, Theor. Chem. Acta 97, 14 (1997). 

W. J. Hunt and W. A. Goddard III, Chem. Phys. Lett. 

3, 414 (1969). 

D. Healion, H. Wang, and S. Mukamel, 

|J. Chem. Phys . 134, 124101 (2011)' 

J. D. Biggs, Y. Zhang, D. Healion, and S. Mukamel, J. 

Chem. Phys. 136, 174117 (2012). 

J. H. D. Eland, M. Tashiro, P. Linusson, M. Ehara, 
K. Ueda, and R. Feifel, Phys. Rev. Lett. 105, 213005 
(2010). 

N. V. Kryzhevoi, R. Santra, and L. S. Cederbaum, J. 
Chem. Phys. 135, 084302 (2011). 

L. Inhester, C. F. Burmeister, G. Groenhof, and 
H. Grubmiiller, J. Chem. Phys. 136, 144304 (2012). 
R. Santra, N. V. Kryzhevoi, and L. S. Cederbaum, 
Phys. Rev. Lett. 103, 013002 (2009). 



13 



[108] J. Schirmer, L. S. Cederbaum, and O. Walter, Phys. [124 
Rev. A 28, 1237 (1983). 

J. Schirmer and A. Barth, Z. Phys. A 317, 267 (1984). [125; 
J. Niskanen, P. Norman, H. Aksela, and H. Agren, J. 
Chem. Phys. 135, 054310 (2011). [126; 
O. Takahashi, M. Tashiro, M. Ehara, K. Yamasaki, and 
K. Ueda, J. Phys. Chem. A 115, 12070 (2011). [127; 
O. Takahashi, K. Yamasaki, S. Nagaoka, and K. Ueda, 
Chem. Phys. Lett. 518, 44 (2011). [128 

A. T. B. Gilbert, N. A. Besley, and P. M. W. Gill, J. [129^ 
Phys. Chem. A 112, 13164 (2008). 

U. Ekstrom, P. Norman, V. Carravetta, and H. Agren, [130 
Phys. Rev. Lett. 97, 143001 (2006). 

M. Linares, S. Stafstrom, Z. Rinkevicius, H. Agren, and [131 
P. Norman, J. Phys. Chem. B 115, 5096 (2011). 

B. Brena, P. E. M. Siegbahn, and H. Agren, J. Am. [132 
Chem. Soc. 134, 17157 (2012). [133^ 
A. Nardelli, G. Fronzoni, and M. Stener, Phys. Chem. 
Chem. Phys. 13, 480 (2011). [134 
J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 [135^ 
(1981). 

G. Tu, Z. Rinkevicius, O. Vahtras, H. Agren, U. Ek- 
strom, P. Norman, and V. Carravetta, Phys. Rev. A 
76, 022506 (2007). [136; 
A. Nakata, Y. Imamura, T. Ostuka, and H. Nakai, J. 
Chem. Phys. 124, 194105 (2006). [137; 
A. Nakata, Y. Imamura, and H. Nakai, J. Chem. Phys. 
125, 064109 (2006). [138; 
Y. Imamura, R. Kobayashi, and H. Nakai, J. Chem. 
Phys. 134, 124113 (2011). [139; 
J.-W. Song, M. A. Watson, A. Nakata, and K. Hirao, 
J. Chem. Phys. 129, 184113 (2008). [140; 



109 
110^ 

111 

112 

113; 

114 

115 

116; 

117; 

118; 

119 



120 



121 



122 



123 



D. Prendergast and G. Galli, Phys. Rev. Lett. 96, 
215502 (2006). 

M. Taillefumier, D. Cabaret, A.-M. Flank, and 

F. Mauri, Phys. Rev. B 66, 195107 (2002). 

A. Ipatov, F. Cordova, L. J. Doriol, and M. E. Casida, 

J. Mol. Struct.: THEOCHEM 914, 60 (2009). 

Z. Li, W. Liu, Y. Zhang, and B. Suo, J. Chem. Phys. 

134, 134101 (2011). 

Z. Li and W. Liu, J. Chem. Phys. 135, 194106 (2011). 
1. V. Schweigert and S. Mukamel, 
|Phys. Rev. A~77, 33802 (2008) , 

J. D. Biggs, J. A. Voll, and S. Mukamel, Phil. Trans. 
R. Soc. A 370, 3709 (2012). 

O. Roslyak, B. Finger hut, K. Bennet, and S. Mukamel, 
(Submitted, 2012). 

A. D. Becke, J. Chem. Phys. 98, 5648 (1993). 

P. Stephens, F. Devlin, C. Chabalowski, and M. Frisch, 

J. Phys. Chem. 98, 11623 (1994). 

T. H. Dunning, J. Chem. Phys. 90, 1007 (1989). 

M. Valiev, E. Bylaska, N. Govind, K. Kowalski, 

T. Straatsma, H. van Dam, D. Wang, J. Nieplocha, 

E. Apra, T. Windus, and W. de Jong, Comput. Phys. 
Commun. 181, 1477 (2010). 

J. Stohr, NEXAFS Spectroscopy (Springer- Verlag, 1992) 
pp. 34-39. 

Z. Li, D. Abramavicius, and S. Mukamel, J. Am. Chem. 
Soc. 130, 3509 (2008). 

S. Mukamel, Principles of Nonlinear Optical Spec- 
troscopy (Oxford University Press, USA, 1999). 
D. L. Andrews and T. Thirunamachandran, 
|J. Che m. Phys. 67, 5026 (1977), 

M. Robin, I. Ishii, R. McLaren, and A. Hitchcock, J. 
Electron. Spectrosc. Relat. Phenom. 47, 53 (1988). 



